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FOREWORD 

The  Summary  Report  was  prepared  by  RCA  Laboratories,  Princeton, 

New  Jersey,  under  Contract  No.  F30602-72-C-0486 .  It  describes  work  performed 
' r°m  June  1972  to  March  1973  in  the  Communications  Research  laboratory, 

Dr.  K.  H.  Powers,  Director.  The  principal  investigator  and  project  scientist 
is  Dr.  D.  A.  de  Wolf. 

The  report  was  submitted  by  the  author  in  April  1973.  Submission 
of  the  report  does  not  constitute  Air  Force  approval  of  the  report's  findings 
or  conclusions.  It  is  submitted  only  for  the  exchange  and  stimulation  of 
ideas . 

The  Air  Force  Program  Monitor  is  Captain  Darryl  P.  Greenwood. 

This  technical  report  has  been  reviewed  and  is  approved. 


SUMMARY  AND  ABSTRACT 


A  summary  Is  given  of  the  results  up  to  February  1973  obtained  under 
Contracts  F30602-71-C-0356  and  F30602 -72-C-0486  regarding  the  behavior 
of  laser  beams  in  turbulent  air.  It  appears  important  to  make  a  distinction 
between  focused  and  other  beams.  Amplitude  scintillations  appear  to  be 
unimportant  in  the  former.  Angle-of -arrival  power  spectra  have  been 
computed  for  a  simple  interferometer  and  these  also  contain  the  temporal 
information  of  focal-spot  fluctuations. 
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L  Lo  LB  Lo 

focal  spot  of  a  Gaussian  beam. 
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1.  INTRODUCTION 


Results  obtained  under  the  contract  efforts  to  RADC,  numbered 
F30602 -7 1-C -0356  and  F30602-72-C-0486 ,  in  a  two-year  period  from  June 
1971  to  August  1973  have  been  reported  at  regular  intervals  in  a  number 
of  reports,  TRI-TRVU  [  1  ]  .  These  results  are  therefore  somewhat  dispersed, 
and  -  moreover  -  several  of  them  have  required  modification  in  later  re¬ 
ports  from  earlier  versions.  This  special  report  serves  as  a  summary  of 
most  of  the  results  obtained  up  to  a  time  in  early  1973,  and  it  includes 
some  new  results  yet  to  be  reported  in  more  detail. 

"he  purpose  of  the  contract  effort  is  to  analyze  the  degradation 
of  laser  beams  in  air  by  turbulence.  Clear -air  turbulence  manifests  it¬ 
self  as  irregular,  small  variations  in  temperature  over  a  large  area  of 
the  ground  (perhaps  tens  of  kilometers).  The  refractive  index  is  a  func¬ 
tion  of  temperature,  and  hence  it  fluctuates  irregularly  causing  rays 
in  the  beam  to  deviate  from  straight-line  paths.  Angular  fluctuations 
(beam  wander  spread),  and  scintillations  in  the  image  plane  (blurring, 
distortion)  are  the  result. 

The  approach  of  the  work  on  these  contract  efforts  has  been  to  give 
a  statistical  description  of  these  various  degradations  of  a  laser  beam 
in  terms  of  the  statistics  of  the  underlying  refractive-index  fluctua¬ 
tions  (which  in  turn  are  given  by  temper ature -flue tuat ion  statistics 
that  can  be  measured  in  principle).  Thus,  an  atmospheric  model  is  re¬ 
quired.  We  have  used  the  following  model  for  homogeneous  turbulence  at 
altitude  x.  Let  the  dielectric  permittivity  deviate  from  its  average 
(very  close  to  unity)  by  an  amount  5£  (r'.t).  Define  a  Fourier  transform 
of  the  covariance  of  the  dielectric  permittivity! 


•  1  • 


(1) 


*00 


5/d3ir  (s£(r,t)56(r  +  Ar  ,t)^  exp(iK»Ar) 


where  r  is  a  point  at  altitude  x.  The  atmospheric-turbulence  model  spec- 
cifies  a  spectrum, 


0(K)  =  32r3  X  0.033  C  2(x)(K2+L  "2)"ll/6exp(-K2/K  2)  , 

no  m 


(2) 


where  macroscale  L  is  presumably  a  function  of  x,  microscale  2  is  re- 
o  o 

-1  2 

lated  to  K  by  the  relationship  2  =  5.92 K  ,  and  C  (x)  is  the  refractive 

m  o  m  n 

index  structure  constant  related  to  its  sea-level  value  [2]  by  the  rela- 
t ionshin, 


C  2(x)  =  C  2(x  ) (T  /T)2(x  /x)2/3(1  +  x/2  )"2/3 
n  n  o  o  o  s 

X  exp[ -2(x-xQ)/h]  , 


(3) 


where  xq  «  lm  and  are  sea-level  values  of  effective  altitude  and  tempera 

ture,  and  -7£  is  the  Monin -Obukhov  length  (2  <»  at-  dawn  and  dusk,  and 

s  s 

2  1.5m  at  midday  for  well -developed  turbulence).  The  structure  func- 

s 

2 

tion  (z)  can  be  obtained  from  temperature -variation  measurements, 

and  the  microscale  K  3,  although  harder  to  estimate,  enters  into  the 

m 

calculations  in  such  a  manner  that  an  error  of  a  factor  twr  r  so  in  its 

magnitude  is  not  serious.  The  macroscale  L  varies  more  but  it  often 

o 

drops  out  of  the  end  result. 


A  useful  form  for  the  electric  field  of  a  laser  beam  is  [3], 

E(?)  =  L  3exp[ik(L  +  p2/2L)] 

2  (4) 
X  P1 B( *1  )Uo Cr\ )exp[~Y“  <£  -  |)  -tkp-^/L]  . 

Here  r^  =  (p^,0)  is  a  point  in  the  plane  of  the  transmitting  aperture, 


-  2  - 


and  r  =  (p,L)  is  in  the  plane  of  the  receiver;  U  (r  )  is  an  aperture 

2  2 

or  pupil  function  [e.g.,  exp(-p^  /2rQ  )  for  a  Gaussian  beam];  R  is  the 
radius  of  curvature  of  the  phase  front  at  z  =  0  (R  >  0  for  a  focusing 
beam),  and  B(r,r1>  is  the  ratio  of  the  electric  field  in  turbulent  air 

-X 

at  r  due  to  a  point  source  at  r  to  a  similar  field  in  free  space.  This 

expression  requires  kL  »  1  (radiation-zone  condition),  L  «  kr  2  (near 

field  of  the  iperture),  and  a  sagittal  condition  kr  ^  «  IT  ,  all  of 

which  usually  hold,  and  it  is  restricted  also  to  points  r  such  that 

43 

kp  «  L  which  is  sufficiently  liberal  to  include  all  points  of  interest 
in  the  image  plane  when  propagation  distances  exceed  a  few  hundred  meters. 

Our  work  shows  that  the  effect  of  atmospheric  turbulence  upon  laser 
beams  is  importantly  different  in  the  focal  plane  of  an  optimally  focused 
beam  from  the  effect  upon  coll imated-beam  and  diverging-beam  images.  It 
makes  sense  to  distinguish  the  case  of  a  focused  beam  from  the  others 
because  there  appear  to  be  onlv  negligible  amplitude  effects  in  the  former. 
A  strong  case  can  be  made  that  the  focal  image  of  a  focused  beam  is  de¬ 
termined  by  free-space  diffraction  effects  upon  which  ray  bending  due 
to  turbulence  is  superimposed.  Ensuing  results  are  discussed  in  the 
next  section. 


2  RESULTS  IN  THE  FOCAL  PLANE  OF  A  FOCUSED  BEAM 
2.1.  Log  -Ampl itude  Fluctuations 

The  ele^cric  field  in  the  focal  plane  of  a  focused  beam  is  given 
by  Eq .  (4)  in  the  special  case  R  =  L.  For  a  Gaussian  beam  in  free  space, 
Eq .  (4)  yields  [3], 


(5) 


E  (r)  =  "Lk'  ■  exp(ik[L  +  p2/2(L-ikw2) ] } 

L-ikw 

-2  -2 

w  =  r  4-  ik/R 
o 

Using  methods  of  our  own  that  lead  to  an  intermediate  result  equivalent 
to  Eqs .  (28)  and  (29)  of  Ishimaru  [4]  we  find  a  log-amplitude  variance 
for  a  Gaussian  beam: 

(Xgb2)  *  0.102  (L/krQ2)7/6  X  (Cn2k7/6LU/6)  (fa) 

The  constant  in  Eq  .  (fa)  is  (0  .033it2/12 )  X  T(3  )T(2/3 )/ T(  1 1/3)  ~  0.102. 

2 

The  important  point  is  that  L  «  kr^  for  all  cases  of  interest  (whereas 
0.31C  2k7^L*^6  appears  never  to  exceed  a  value  ot  the  order  of  10), 
and  therefore  <  Xg^  y  «  1,  Thus,  amplitude  fluctuations  do  not  contrib¬ 
ute  appreciably  to  the  effects  in  the  focal  plane. 


2.2.  Focal -spot  broadening 

The  previous  result  strongly  suggests  that  the  image  in  the  focal 
spot  is  distorted  by  undulations  in  the  rays  directed  by  the  mirror  to¬ 
ward  a  focal  spot.  There  are  many  different  ways  to  define  variables 
for  the  broadening.  Some  of  these  are  discussed  in  TRI  and  TRIV .  In  our 

view,  the  most  fundamental  (although  not  necessarily  the  most  practical) 

2 

one  is  the  variance  width  r  .  Consider  the  normalized  average  inten¬ 
sity  i(p)  =  <I(o, L)  /fd2p  I(p,L)>  .  The  variance  width  r^2  is  then 
defined  as 


rL2  =  yd2pp2i(p)  -  [Jd2ppiCo)]2 


(7) 


The  quantity  r  is  defined  likewise  in  terms  of  the  normalized  free- 
space  intensity  iQ(p)  equivalent  in  definition  to  iCp)  if  <I(fs,L)^  is 


-  h  - 


2  2  2 

replaced  by  I  (p,L).  We  finally  define  r  h  r  -rT  ,  and  find: 
°  L 15  L  Lo 


r  2  =  12rt2  X  0.033L2k 
LB  m 


inf ds 


(s) 


(8) 


for  an  axial  path  from  s  =  0  to  s  =  I  .  For  a  homogeneous  path  one  finds 
2  3  1/3  2 

rT B  ~  1.3L  K  Cn  .  The  significance  of  this  result,  which  may  be 

somewhat  difficult  to  compare  in  detail  to  data  because  it  requires  much 

processing  of  data,  is  that  it  is  a  very  fundamental  and  general  result, 

independent  of  beam  shape,  aperture  function,  etc.,  and  it  holds  even 
2  2 

when  r^  and  r^Q  are  not  finite  in  which  case  it  can  be  obtained  from 
data  by  forming, 


2  . 

rLB  J 


d2pp2[i(p)-io(p)] 


(9) 


where  p^  is  a  radius  value  beyond  which  the  intensity  drops  below  the 

2 

noise.  The  quantity  r  represents  an  area  that  expresses  a  displace¬ 
ment  of  energy  away  from  the  central  axis  into  the  sides  by  turbulence- 
induced  refraction.  This  view  is  supported  by  looking  at  the  radial 
displacement  r  '  =  LO  where  0  is  the  angular  deviation  of  a  ray  pointed 
initially  at  the  focal  point  from  some  location  on  the  mirror  and  de¬ 
flected  randomly  and  slowly  by  re fractive -index  fluctuations.  It  is 
easily  seen  that  (o)  =  0  and  <T  02 )  =  I^LkJ^C^2  for  a  homogeneously 

turbulent  medium,  as  a  consequence  of  which  it  follows  that 
2  2 

( ( rLB * )  )  =  rLB  •  This  tesult,  coupled  with  Eq .  (6),  suggests  a  sim¬ 

plified  model  for  the  instantaneous  intensity  in  the  focal  plane  that 
will  be  considered  in  future  work. 

Equation  (8)  has  been  worked  out  for  air-to-ground  propagation  as 


well  as  for  the  opposite  case  (unfortunately,  the  results  in  TRIV  have 


been  interchanged  for  those  two  cases).  In  both  cases,  we  name  the  re  - 

2 

suiting  variance  width  r^  (Z)  where  Z  is  the  height  of  the  receiver 

(ground-to-air  case)  or  of  the  transmitter  (air-to-ground  case).  It  is 
2 

compared  to  r,  (0):  the  horizontal -propagat ion  ^asc.  The  following 

2  2 

table  summarizes  results  for  the  ratio  r  (Z)/rT  (0). 

LB  LB 


TABLE  I : 


— 

Sunny  Day 

1 - — 

Dawn ,  Du  sk 

Ground 

-To-Air 

9x  1/3B  *'Vl 

o  a 

9(x  /Z)2/3 
o 

1 

t 

Air -To 
-Ground 

l 

1.80c  f  )2/V4/3  : 

o  s 

1.3(x  /Z)2/3 
o 

These  formulas  arc  approximations  for  altitude  Z  much  less  than  h 

(see  Eq  .  v3)  for  explanation  of  *  and  x  ).  The  inverse  ratio  has  been 

s  o 

plotted  in  Figs.  1  and  2  for  appropriate  choices  of  £  (1.5  mor  °o)  and 

s 

XQ(lm).  The  calculations  h*ve  been  performed  by  inserting  Eq .  (3)  into 

Eq .  (8),  and  they  include  some  corrections  of  the  original  presentation. 

2  2 

The  ground-level  variance  width  r^  (0)  =  rLfi  has  bnen  plotted  in  Fig. 

2 

3  for  several  values  of  .  Fig.  3  also  includes  two  plots  of 
2  2 

r^Q  =  (L/kr^)  ,  the  variance  width  of  the  free-space  focal  spot  of  a 
diffraction-limited  Gaussian  beam, for  comparison.  It  indicates  the  sig¬ 
nificance  of  the  broadening. 

2.3.  Focal -spot  Scintillation  Rates:  Power  Spectrum 

2 

fhe  variance  width  r^^  ,  defined  in  (7),  can  be  made  into  a  covariance 
finctionwith  respect  to  time  (thus  to  yield  info;nnation  on  scintillation 


-  6  - 


SLANT  PATH  :  UP 


Figure  1.  Reduction  factors  for  the  average  focal-spot  area  in  ground' 
to-air  propagation.  Receiver  at  altitude  Z. 


rates)  by  redefining 


rLB2(T)  =/d2^P2[<E(r,t)E*(?,t+T)>  -  Io(^]//d2p<I(?)>  . 


(10) 


Note  that  the  first  moments  are  ignored  (asymnetric  beams  are  not  con¬ 
sidered).  For  a  homogeneously  turbulent  medium,  it  can  be  seen  that 


rLB2(r)  -  L2  <0(t)G(t+'i)  > 


(11) 


where  fl(t)  is  the  angular  deviation  of  a  ray.  Thus,  the  Fourier  trans- 

2 

form  of  rLB  (t)  yields  the  power  spectrum  of  either  the  square  root  of 
the  focal-spot  variance  width,  or  -  equivalently  -  of  the  transverse 
deviation  of  a  ray.  The  power  spectrum  W  (o))cL  tells  you  how  this  vari¬ 
ance  is  apportioned  into  contribution!  with  fluctuation  rates  lying  be¬ 
tween  (0  and  -o  +  do.  It  has  been  computed  for  turbulence  characterized 
by  ^  random  Gaussian  distribution  of  velocities  with  mean  U  and  variance 
AU  .  Only  the  transverse  components  UT  of  U  and  2AU/3  enter  into  the 
calculations.  The  results  are  lengthy  and  some  preparation  is  needed 
to  write,  them  down.  The  velocities  l'T ,  AU,  and  the  scales  L  ,  SI  define 
four  frequencies 


o  o 


oif  -  Ut/Lc 


=  KmUT 


An  =  4AU/L  V~ 3 
o 

A1  =  4K  ADA/" 3 
m 


(12) 


For  the  turbulence  spi  trum  of  Eq .  (2)  we  find  the  following  power  spec- 
trum: 


W0(oJ)  =  4tt2  X  0.033Cn2LLo"1/3[v}(o),5/6)  -  WM (oo,  11/6)]  (13) 

The  forms  for  W  (<o,q)  can  be  found  from  Tables  II  and  III,  and  they  can 
be  adopted  to  Eq .  (13)  by  setting  q  =  1 l/b  and  subtracting  the  result- 
ing  form  from  that  for  q  =  5/6.  A  normalized  form  of  the  spectrum  based 

10  . 
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TABLE  II:  ANGLE -OF -ARRIVAL  POWER  SPECTRUM 
FOR  FROZEN  FLOW  (£U=0) 


W  U ,  q ) 


R  2.  2  \  l/2-q  -J/nl 


IT  i+u  /ajT 


i  U  [  1/2  ,  3/2-q,  (co2-hu2) /ft2)  ] 


<£  /  2  2  V/2-«  '  ^  i  (l^2/^)'^'"2^ 

r(q)  ^  \1+“  /l“t)  e  T  wx  WT  1 


I 


r (q— 1/ 2)  /£  r (q-l/2  /n 

r  (q)  0JT  r(q'.  «JT 


WT  =  u  /L 

I  To 


(«  L  )  -  P- 

mo  o)T  ^  wT 


=  KmUT 

1  m  T 


9  9 

-2q  -u)“/i2T 
e 


WA  (c  ,q) 


TABLE  III:  ANCLE -OF -ARRIVAL  POWER  SPECTRUM 
FOR  RANDOM  FLOW  (UT=0) 


r  (q-l/2)  U(q-l/2,  1/2  qj2/ A  to 2 ) 


^  O,  2  1/2_q 

Ac  2(KmL0w/Aa,)  Kq-i/2  (2w/A^ 


—(^.1/2)  JL  r  (  1/2)(  H  \  ^  ,  _rr_  / 2w  \ 

(q)  Aui  \  Aw  j  Ao  K  mL0^  Aw  (  ASI  j 


-q  -2io /Aft 


Au  -v  AU/L 


Aft  9,  <  AU 


IPPI 


on  Table  II  is  plotted  in  Fig.  4  for  K  L  =  100  and  K  L  =  1000.  All 

mo  mo 

frequencies  denoted  by  letter  "f"  in  this  plot  can  also  be  considered 
as  angular  frequencies  U3. 

The  power  spectrum  has  a  flat  portion  for  low  frequencies,  followed 

-2/3 

by  a  transition  to  an  03  portion,  which  in  turn  is  followed  by  a 

-2/3 

rapid  decay.  The  03  portion  is  presumably  the  only  significant  part 
of  the  spectrum  because  the  lower-frequency  portion  assumes  the  absence 
of  large-scale  fluctuations,  and  the  high-frequency  fall -off  is  depen¬ 
dent  in  detail  upon  the  exponential  tail  of  Eq .  (2):  an  approximate 
description  which  suffices  for  most  other  applications  but  not  for  this 
one.  The  following  approximation  of  W  (to, 5/6)  appears  to  yield  a 
useful  expression  whenever  0)  exceeds  the  maximum  of  cqj.  and 

W  (03,5/6)  «  r(l/3)  ^  (S7)"2/31F1(1/6,1,  -o^W)  ,  (14) 

where  ^F^  is  a  confluent  hypergeometric  function  in  a  well-known  nota¬ 
tion  [7].  An  asymptotic  expansion  for  small  Ato/o^  is  probably  sufficient 
for  all  cases  except  those  where  UT  «  0 . 

The  spectrum  in  Fig.  4  has  been  tested  with  Russian  data  [6]  shown 
in  a  thick  interrupted  line.  These  are  phase-difference  data,  not  angle- 
of-arrival  data,  and  therefore  agreement  is  good  only  for  f  <  UT/d  where 
d  is  the  separation.  For  these  data,  that  implies  f/ f^,  <  10. 

3 .  RESULTS  FOR  COLLIMATED  (PLANE -WAVE)  AND  DIVERGING  (SPHERICAL-WAVE) 
BEAMS 

3.1.  The  Variance  Width  At  z  =  L 

For  a  Gaussian  beam  described  by  (5),  one  can  recompute  Eq .  (7) 
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f 


to  find: 


rT2  =  r  2 [ (1 -z/R)2  +  (z/kr  2)2]  +  r.  2  ,  (15) 

L  O  O  LB 

i.e.,  the  variance  width  in  any  plane  at  an  appreciable  distance  from 

the  laser  is  the  Pythagorean  sum  of  the  width  r  solely  determined  by 

atmospheric  turbulence  and  a  width  rT  characteristic  of  the  beam  in 

Lo 

free  space.  At  z  =  R,  Eq .  (15)  reduces  to  the  focused-beam  case  in  the 

focal  plane  depicted  in  Fig.  3.  For  a  collimated  beam,  it  reduces  to 

2  2  2  22222 
rT  «  r  “  +  rT  „  ,  and  for  a  diverging  beam  to  rT  «  r  (R  +  Z)  /R  +  rT  „  . 
L  O  LjD  L  O  Lj  15 

Thus  for  wide  beams,  or  for  diverging  beams,  the  broadening  effect  is 
much  less  important  than  it  is  in  the  focal  plane  of  a  focused  beam. 

The  internal  structure  of  the  large  rpot  at  z  =  L  is  now  of  much  greater 
importance,  and  it  appears  to  be  de ten  lined  mainly  by  amplitude  fluctu¬ 
ations:  in  strong  contrast  to  the  focused-beam  case. 

3  .2  .  Log-amplitude  Fluctuations 

The  analysis  of  amplitude  fluctuations  -  probably  the  dominant 
mechanism  of  image  distortion  in  diverging  beams  -  is  extremely  diffi¬ 
cult.  The  following  results  have  been  obtained,  and  they  are  illustrated 
in  Fig.  5  for  plane  waves  (the  illustration  for  spherical  waves  would 
be  quite  similar). 

3.2.1.  Weak -amplitude  regime 

The  amplitude  fluctuations  are  log-normally  distributed  with  zero 
log-amplitude  mean,  and  with 

{  I®)  =  I  N  exp[2N(N-l)<'X2)  ]  (16) 

o  r 

where , 


-  U  - 


LOG-NORMAL  II  RAYLEIGH  I 


S)^C  *  *  /  7/3.  3^2  -1/6  I 
/J^°-4(ACm  ^  Cn) - L - -  — — — 


(k/K~  L )  • 


(kLc0/L) 


Fipure  5.  The  plane-wave  log-amplitude  variance  <x>  as  a  function 
_  2  .  „„  7/6.  ll/6„  2  .  .  . 


0.31k7/6L11/6C  2. 

n 


A  sample  curve  has  been  sketched. 


Note  that  the  two  abscissa  points  really  correspond  to  values 
of  kL^K  /3C  which  does  not  differ  greatly  from  O  ^ . 


:  plane  waves 


07) 


<X2>  ~  0.307  k7/6LU/6Cn2 
<(x2>  «  0.124  k7/6LU/6C  2 


spherical  waves 


These  modif ied-Rytov  results  appear  to  be  valid  whenever  the  variance 

1/3  2 

of  the  ray-angle  deviation  K  LC  is  much  less  than  the  square  of 

m  n 

the  Fresnel  angle  (kL)  J^2  i.e.,  kL2.<  *^3C  2  «  1. 

m  n 


3.2.2.  Saturation  regime 

2  1/3  2 

As  the  parameter  kL  exceeds  unity,  i.e.,  when  ray  bending 

starts  to  lead  to  interference  effects  with  undeflected  rays,  the  vari¬ 
ances  in  Eqs.  (17)  are  no  longer  small  compared  to  unity  and  Eqs.  (17) 
are  no  longer  valid.  This  parameter  regime  is  quite  relevant  to  laser - 
beam  propagation  over  kilometer  paths.  New  results  indicate  that  the 
amplitude  is  asymptotically  log-normal,  i.e.,  Eq .  (lb)  holds  asymptotically. 
Recent  results  by  other  workers  [5]  indicate  that  the  amplitude  cannot 
be  exactly  log-normal;  however,  an  asymptotic  result  seems  resonable 
at  present,  and  Eqs.  (17)  must  then  be  replaced  by, 


waves 


<x2>  -  0.41(/<m7/3L3Cn2)“1/6  :  plane 

:  spherical  waves 


(y2)  *  ^  7/3t3 „  2, -1/6 


(18) 


1.25(k  ,/jLjC  )' 
m  n 


Thus,  the  amplitude  fluctuations  have  peaked  somewhere  at 
1/3  2  -1 

Km  LCn  ~  *  and  now  decrease  with  increasing  L  and  C  2. 


3.2.3.  Rayleigh  regime 

The  ^crease  does  not  continue  ad  infinitum  because  when  kL2K  *^3C  2 

m  n 

approacnes  and  exceeds  the  Fresnel  number  kL  2/L,  the  amplitude  becomes 
Rayleigh  distributed,  i.e., 


<iN>  =  n!i 


N 


(19) 
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This  is  a  limit  never  observed  with  laser  beams,  because  the  above  cri¬ 
terion  means  that  rays  are  refracted  laterally  by  more  than  a  macroscale 

Lo’  and  that  is  very  unlikely  for  beams  propagating  even  tens  of  kilometers 
through  turbulent  air. 


4 .  INTER!”  ROMETER;  FEATURE  RESOLUTION 
4 • 1  •  Resolution  Of  Features  Of  An  IUuminated  Object 

A  simple  study  has  led  to  a  criterion  for  distinguishing  two  fea¬ 
tures  separated  by  an  angle  0  =  p  /L  on  an  object  at  distance  L  from 

o  ^ 

the  observer  in  the  presence  of  turbulent  air.  Complete  absence  of  other 
noise  has  been  assumed  so  that  perfect  resolution  at  any  angle  in  free 
space  i s  assumed.  The  criterion  expresses  the  magnitude  or  the  varia¬ 
tion  50  in  the  difference  of  two  rays  to  the  eye  coming  from  each  point 
feature .  It  reads 

<((5 0/9  )2y&  (r  JO. 6U  )2  for  p  «Z 

g  L.H  O  O  O 

2  (20) 

~  (2  •4rT  J P  )  for  p  »  Z 

Lb  O  O  O 

Thus,  small  separations  compared  to  Z  are  either  obscured  (r  >  £  } 

°  LB  o ' 

or  they  are  not.  Larger  separations  can  be  resolved  or  not,  depending 
upon  the  size  of  rLg  (to  be  read  off  from  Fig.  3  for  horizontal  propa¬ 
gation).  Slant-path  propagation  requires  some  modification  factors, 
namely, 

<  5A2(Z)>  /(bO2)  «  5(zo/Z) 

*  3(z  /Z)2/3 
o 

for  an  object  at  altitude  Z. 

4*2.  Power  Spectrum  Of  Interferometer  Difference  Angle -Of -Arrival 

Finally,  a  simple  application  and  extension  of  the  results  of  sec 
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sunny  day 
dawn,  dusk 


(21) 


tion  2.3  are  to  be  found  in  examining  the  de  ference  angle  of  arrival 
of  a  simple  interferometer  consisting  of  two  parallel  rays  separated  by 
a  distance  p .  This  angle  is  a  measure  of  the  distortion  of  an  image, 
whereas  the  ha^f-sum  angle  expresses  a  lateral  shift  of  the  image.  The 
power  spectrum  for  this  case  has  a  form  similar  to  Eq .  (13),  namely 

2  2-1/3 

W.  (u))  =  8n  X  0.033C  LL  '  [wft  (ca,5/6)  -  Wcn(ca,  11/6)]  .  (22) 

otf  n  o  ot7  5" 

The  difference  lies  in  the  form  of  W  , 

OH 


w5f)(“»q)  =  W0(co,q)F  (ce.q)  ,  (23) 

which  differs  from  the  single-ray  function  W#(ce,q)  by  inclusion  of  a 
filter  factor.  Some  extra  definitions  are  required  in  order  to  give 
the  filter  factor.  Let  cp  be  the  angle  in  the  x-y  plane  between  trans- 
verse  average  velocity  UT  and  separation  vector  p.  Define 

5  -  pcosC?/L 

•  IT0  (24) 

I  =  psin'p/L 
s  o 

The  following  forms  have  been  found  for  the  filter  factor  F  (co,q): 


4.2.1.  Frozen  Flow  (AU  =  0) 

00 

F0(a',q)  =  1  -  £cos(W  ~t  [-(I  2/4)(l-ho2/co  2)]m 

m=u 

U[l/2  +  m,  3/2-q+m,  (ca2-*^2  ] 

UEl/2,  3/2 -q,  (cu2-^2)/^2] 


(25) 


where  U(a,b,z)  is  the  Krummer  function  of  formula  (13.1.33)  of  Ref.  [7]. 
If  the  dissipation  range  of  turbulence  is  excluded,  i.e.,  ca  «  fl  Eq  . 
(25)  is  considerably  simplified  to: 
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(2b) 


23/2-q 


Fp(aMl)  =  1  ■  iwTTl)  cos^u/^Ul  +a?/^2f/2-l/u^-i-1/2 

XKl/2-q[^s(1  +“2/'^)1/2] 


4.2.2.  The  g>2^  portion 

Analogous  to  Eq .  (14),  we  are  able  to  work  out  a  general  expression 
for  the  inertial  subrange  of  the  power  spectrum,  that  is  most  useful 
when  MT  «  UT .  Only  the  q  =  5/6  fi.ter  factor  will  be  given  to  corres¬ 
pond  with  Eq  .  (14)  according  to  the  recipe  of  Eq .  (22):  the  q  =  11/6 
factor  can  be  given  but  it  appears  superfluous  to  do  so: 


2cos|  a>/cO  °° 

(-.5/6>  -  1  - . T(1/3)  £  ir  «  0./2V  /3 

m=o 

iF  (l/6+m,l,-(n  ' 2 /to2'1 

*  W  W  T(7/6  ,  h  h 

1F1(1/6,1  ,-0^  /co  ) 


The  flat  portion  of  the  spectrum  (for  very  small  to)  can  also  be  given  - 
see  TRVI,  section  * .2 .2  -  but  it  is  not  worth  repeating  here  because  it 
is  probably  only  a  formal  result  without  much  physical  significance. 


^•2.3.  Effect  ot  i no  filter  factor 

Equations  (25) -(27 )  are,  perhaps  unnecessarily,  complicated,  but 
their  content  can  be  stated  more  simply  in  ruaHtative  terms.  Just  con¬ 
sider  the  cases  where  F  (o),q)  reduces  either  to  F  1  or  F  -  0  Tn 

P  o 

the  case  that  0,  i.e.,  that  the  plane  of  the  interferometer  is  per¬ 

pendicular  to  the  crosswind,  it  follows  f rom  Eq .  (26)  that  F  ~  1  for 

2/1  P 

frequencies  well  into  the  a)  region.  Thus,  for  a  vertical  interfer¬ 
ometer,  the  filter  factor  is  important  only  at  low  frequencies.  On  the 
other  hand,  for  a  horizontal  interferometer,  i.e.,  £  —  0,  the  filter 

factor  reduces  to  l-cos^oi/u^)  for  all  except  the  very  highest  frequencies 
(depending  upon  how  non-zero  really  is).  The  filter  factor  < h.  ,  filters 
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out  of  the  spectrum  the  resonance  frequencies  u>  =  2jtna4j,/?c  =  2jrnU^,/pcosCD 
~  2"nU^,/p,  but  does  not  otherwise  affect  the  spectrum. 

5.  CONCLUDING  REMARKS 

The  preceding  sections  2-4  summarize  the  most  important  results  of 
the  contract  effort  to  date.  Ir.  summary,  it  appears  that  a  focused  beam 
produces  an  image  spot  in  turbulent  air  that  is  determined  by  free-space 
diffraction  interference  and  turbulence -induced  angular  refraction  effects. 
This  result  is  expressed  by  Eq .  (6)  for  negligible  turbulence -induced 
amplitude  fluctuations,  and  by  Eq .  (8)  for  the  average  area  broadening 
of  the  focal  spot.  The  latter  broadening  corresponds  exactly  to  a 
broadening  due  to  angular  refraction  of  rays  only.  Table  I  expands  Eq . 

(8)  to  the  case  of  slant  paths.  Thus  the  quantity  determining  the  broaden¬ 
ing  of  the  focal  spot  is  in  essence  the  lateral  deviation  of  a  randomlv 
refracted  ray.  The  power  spectrum  associated  with  both  of  these  is  one 
and  the  same  quantity  and  it  is  given  by  Eqs.  (13),  (14),  and  Tables 
II  and  III.  The  power  spectrum  subdivides  a  variance  (in  this  case  of 
either  the  focal  radius  or  the  lateral  deviation  of  a  ray)  into 
frequency  components;  each  component  corresponding  to  fluctuations  at 
a  given  frequency.  It  contains  all  the  information  on  fluctuation  rates. 

Strongly  diverging  or  wide  collimated  beams,  on  the  other  hand, 
produce  images  in  which  amplitude  flucutations  are  also  of  importance, 
particularly  in  the  interior  of  the  spot.  This  information  is  not  con¬ 
tained  in  the  mutual-coherence  factor,  apparently,  because  the  mef  is 
determined  solely  by  phase  fluctuations.  Equations  (16)-(19)  give  re¬ 
sults  for  amplitude  fluctuations.  Not  included  in  these  results  are 


20  - 


the  effects  of  intermittency :  these  will  be  reported  at  a  later  date. 

Some  simple  criteria  for  resolution  of  features  of  an  image  in  tur¬ 
bulent  (but  otherwise  noise-free)  air  are  given  in  Eqs.  (20)-(21).  As 
expected,  these  criteria  revolve  around  angular-deviation  parameters. 

Finally,  the  pcwer  spectrum  of  a  simple  interferometer  difference  - 
angle  of  arrival  is  computed  in  Eqs.  (25),  (26),  and  (27),  and  related 
to  the  angle-of-arrival  power  spectrum.  These  results  complement  those 
of  Reinhardt  and  Collins  [8]  for  phase-difference  spectra. 

The  results  appear  to  contain  much  information.  A  power  spectrum, 
for  instance,  probably  contains  much  more  information  than  one  really 
requires  for  understanding  the  rates  at  which  fluctuations  occur.  One 
of  the  tasks  of  the  immediate  future  will  be  to  extract  pertinent  in¬ 
formation  from  these  many  results. 
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